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ABSTRACT 

We present the source catalog of a new Chandra ACIS-I observation of NGC 300 obtained as part of 
the Chandra Local Volume Survey. Our 63 ks exposure covers ~ 88% of the 1^25 isophote {R ~ 6.3 kpc) 
and yields a catalog of 95 X-ray point sources detected at high significance to a limiting unabsorbed 
0.35-8 keV luminosity of ~ 10^^ erg s~^. Sources were cross-correlated with a previous XMM-Newton 
catalog, and we find 75 "X-ray transient candidate" sources that were detected by one observatory, but 
not the other. We derive an X-ray scale length of 1.7±0.2 kpc and a recent star formation rate of 0.12 
Mq yr~^, in excellent agreement with optical observations. Deep, multi-color imaging from the Hubble 
Space Telescope^ covering ^ 32% of our Chandra field, was used to search for optical counterparts to 
the X-ray sources, and we have developed a new source classification scheme to determine which 
sources are likely X-ray binaries, supernova remnants, and background AGN candidates. Finally, 
we present the X-ray luminosity functions (XLFs) at different X-ray energies, and we find the total 
NGC 300 X-ray point source population to be consistent with other late-type galaxies hosting young 
stellar populations (< 50 Myr). We find XLF of sources associated with older stellar populations has 
a steeper slope than the XLF of X-ray sources coinciding with young stellar populations, consistent 
with theoretical predictions. 

Subject headings: galaxies: individual (NGC 300) — galaxies: spiral — X-rays: galaxies: luminosity 
function 



1. INTRODUCTION 

The X-ray emission from normal, non-active galaxies 
originates from a mixture of stellar sources and hot, dif- 
fuse gas. In spiral galaxies, this X-ray emission is domi- 
nated by X-ray binaries (XRBs) comprising either a neu- 
tron star (NS) or a black hole (BH) accreting from a stel- 
lar companion. The fast evolutionary timescale of mas- 
sive stars (~ 10^ years) makes X-ray emission from high- 
mass XRBs (HMXBs) nearly simultane ous with their for- 
mation (Shtykovskiy & Gilfanov 2007), while the longer- 
lived low-mass XRB (LMX B) systems trace older un- 



derlying stellar populations ([ Kong et al. 2002 Soria 
' Kong||2002[ [Trudolyubov et al.|[2UU^ thusTtlie age and 



star formation history (SI^'H) of the host galaxy should 
correlate with the shape of t he X-ray luminosity function 
(XLF) (Grimm et al. 2003] [ Belczynski et al .||2004l |Era-| 
[cleous et al. 2006). Dependencies of the XKB population 
properties on host galaxy metallicity may provide useful 
constraints on binary evolution and population synthesis 
models. 

With its excellent angular resolution (~0'.'5) and posi- 
tional accuracy, the Chandra X-ray Observatory is the 
only X-ray telescope capable of separating the X-ray 
point source populations of nearby galaxies from dif- 
fuse emission. When combined with deep, optical Hubble 
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Space Telescope (HST) imaging, reliable source identifi- 
cation and optical counterpart identification may be car- 
ried out even at distances of a few Mpc. The Chandra 
Local Volume Survey (CLVS, P.I. Benjamin Williams) 
is a deep, volume-limited X-ray survey of five nearby 
galaxies (NGC 55, NGC 300, NGC 404, NGC 2403, and 
NGC 4214) with matched iJST observations down to My 
-^0 (Dalcanton et al.^2009} . When c ombined with the al- 



Tiillmann et al.L 



of the stellar mass and 
tion out to a distance of 



ready well-studied Jisks of M 3 1 (iKong fc Di Stefano 
20031 and M 33 ([Williams et al. || 20U8[ [Flucinsky et al 

these galaxies contain ^997o 
% of the recent star forma- 
3.3 Mpc (Tikhonov et al.|2005l 



Mafz-Apellaniz et al.||2002t [Freedma n & Madore 19T 
Additionally, these galaxies span a representative saniple 
of disk galaxies with a range of masses, metallicities, and 
morphologies. 

In this paper, we present the results of a new Chandra 
observation of the nearby SA(s)d spiral gala xy NGC 300 , 
located at a distance of 2.0 Mpc (Dalcanton et al.||2009 ) 
in the Sculptor Group. NGC 300 is a near-optical twm of 
the Local Group galaxy M 33, and has a low foreground 



Galactic absorption of 4.09x10^ cm ^ (Kalberla et al. 



2005) and is viewed at an inclination angle i = 42""; Ta- 
ble[l] lists some basic properties of NGC 300, as well as 
M 33 and the Milky Way for direct comparison. 

The first X-ray population studies of NGC 300 were 
performed with ROSAT between 1991 and 1997 ( [Read 
"& Pietsch"'2001), identifying 29 bright X-ray sources m- 
cluding NGC 300 X-1, a highly variable supersoft source, 
and other bright sources coincident with known super- 
nova remnants (SNRs) and H II regions. More recently, 
an XMM-Newton survey of the the X-ray source popula- 
tion of NGC 300 was presented by Carpano et al.| (2005 ) 
down to a limiting luminosity of ^ 3 x lO'^^ erg s~ in the 
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TABLE 1 
Summary of NGC 300 Properties 



Property 


NGC 300 


M 33 


Milky Way 


(1) 


(2) 


(3) 


(4) 


Distance 


2.0 Mpc" 


800 kpc^ 


8.5 kpc 


Type 


SA(s)d^ 


SA(s)cd^ 


SBc 


Mb 


-17.66^ 


-18.4^ 


-20.3^ 


Scale length (kpc) 


1.3^ 


1.4^ 


3.00±0.22^ 


Circular velocity 


97 km s-^f 


130 km s-^g 


239±5 km s-^^ 


Estimated stellar mass 


4.3 X 10^ Mq f 


4.5 X 10^ Mq 3 


6.43 ± 0.63 X 10^0 Mq ^ 



Note. — "IPalcanton et al.| ([2009|>; ^illiams et ak] (|2009t; ^N ED; ^Vila-Costas Edmunds] ( |l992| ; ^Munoz-Mateos et al.] ( |2007t ; 
■/ [Puclie et al.||l990| >; ^Corbelli Salucci] ( |2000t ; ^McMillan| ( |2011| > 

NGC 300, with numerous X-ray sources visible. Soft 
X-ray emission (0.35-1 keV) is shown in red, medium 
X-ray emission (1-2 keV) is shown in green, and hard 
X-ray emission (2-8 keV) is in blue. The 1^25 isophote is 
superimposed in white. 

Nine HST fields were used to search for optical coun- 
terparts of each of our X-ray point sources. Three 
HST Advanced Camera for Surveys (ACS 



0.3-6 keV band. A total of 163 sources were detected in 
the energy range of 0.3-6 keV, and the 86 sources falling 
within the 1^25 optical disk were further characterized 
using hardness ratios. X-ray fiuxes, and the available 
ground-based optical imaging. A brief summary of pre- 
vious X-ray NGC 300 studies is provided in Table [2] 

Before the publication of the full Chandra X-ray point 
source catalog, our collaboration employed this new ob- 
servation of NGC 300 to study two individual objects 
of interest: the "supernova impostor" SN 2010da, which 
we argued is consistent with a Be-X-ray binary (BeXB, 
comprised of a young, rapidly rotating B- type star with 
strong emission lines and a NS) origin (Binder et al. 
201 Ib j), and the BH + Wolf-Rayet binary NGC 300 X- 
1 ("Binder et al. 2011a). These sources are included as 
part of the X-ray source catalog, but are not discussed 
in detail in this work. 

The organization of this paper is as follows: in Sec- 
tion [2] we provide the details of our new observation, 
data reduction procedure, and the technique used to gen- 
erate the X-ray source catalog. An analysis of the X-ray 
source catalog, including comparisons to previous stud- 
ies, hardness ratios, and the radial source distribution 
is given in Section |3] In Section [4j we present a new 
X-ray /optical source classification scheme, identify can- 
didate optical counterparts to our X-ray sources, and as- 
sign likely source classifications (HMXB, LMXB, SNR, 
or background AGN) to each source. The X-ray luminos- 
ity function, and the possible role of transient sources in 
determining the luminosity function shape, is discussed 
m Section^ The paper concludes with a summary and 
the main conclusions in Section [6l 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Observations 

NGC 300 was observed on 2010 September 25 for 63 
ks using ACIS-I during the Chandra X-ray Observatory 
Cycle 12, observation ID 12238. Data reduction was car- 
ried out with CIAO v4.3 and CALDB v4.4.2 using stan- 
dard reduction procedures. Images were created with bin 
sizes of 1, 2, 3, and 4 pixels for each of the following en- 
ergy bands: 0.35-8, 0.35-1, 1-2, and 2-8 keV. Exposure 
maps were constructed using mkinstmap and mkexpmap 
in CIAO. Exposure-corrected images were generated us- 
ing our exposure maps and instrument maps. The in- 
strument maps were created assuming spectral weights 
appropriate for a power law spectrum typical of XRBs 
and background AGN (F = 1.9). The average foreground 
column density A^h was used. 

Figure [l] shows an RGB-rendered X-ray image of 



Ford et al. 

1998 ) observations of NGC 300 were taken as part of 
tEilicS Nearby Galaxy Survey Treasury (ANGST, GO- 
10915) forming a radial strip from the center of the 
galaxy into the disk, hereafter referred to as 'WIDEl- 
3.' Six additional archival ACS images scattered across 
the NGC 300 disk (GO-9492) are also included, and la- 
beled as 'NGC300-1-6.' Detailed descriptions of the HST 
star catalogs, including photometric measurem ents and 
quality cuts, a re giv en in Williams et al. (2009) and |Dal-| 
canton et al. (2009), while a thorough discussion 01 the 
survey completeness and sta r formation histories (S FH) 
of NGC 300 are presented in'Gogarten et al | (| 2Q10l). 

Figure [2] shows a ground-based Palomar Observatory 
Sky Survey (POSS) image of NGC 300, with the loca- 
tions of the HST fields shown in red and our Chandra 
field superimposed in blue. The locations of the Chandra- 
detected X-ray point sources (see ^ are shown by green 
crosses. 

2.2. Image Alignment 

To classify optical counterparts to our Chandra X-ray 
sources, we place both the Chandra and i7ST frames onto 
the International Celestial Reference System (ICRS) by 
finding matches between stars or background galaxies in 
the Two Micron All S ky Survey (2MASS) Point Source 
Catalog (Skrutskie et al. 2006). We performed relative 
astrometry by individually aligning each HST diud Chan- 
dra field to the same reference image frame and coordi- 
nate system. 

For our reference image, we retrieved a publicly- 



available grou nd-based R band image (see Larsen & 
Richtler| [2QQQl for more details) of NGC 300 from NED. 
Foreground zMASS stars with ground-based i?-band 
counterparts were identified, and we used the IRAF task 
ccmap to compute the plate solution and update the im- 
age header with corrected WCS information. By using 
the 2MASS stars as reference, the astrometric solution 
yielded root-mean-square (r.m.s.) residuals of 0^00773 in 
right ascension and 0^0225 in declination relative to the 
2MASS astrometry for the ground-based i?-band image. 

We were next able to align each of our HST frames 
and our Chandra frame to the 2MASS-aligned ground- 
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TABLE 2 

Summary of Prior X-ray Surveys of NGC 300 



Observatory 


Number of Detected 


Limiting Unabsorbed Luminosity 


References 




Sources (within D25) 


(1036 erg s-i)" 




(1) 


(2) 


(3) 


(4) 


ROSAT 


47 (29) 


- 4.0 (-7.7) 


[Read & Pietsch' ('2001 1> 


XMM-Newton 


163 (86) 


^ 0.3 (~0.4) 


ICarpano et al. (2o5^ 


Chandra 


95 (77) 


1.5 


this work 



Note. — "The hmiting unabsorbed luminosities correspond to the 0.1-2.4 keV and 0.3-6.0 keV energy ranges for ROSAT and XMM- 
Newton, respectively. In parenthesis we provide expected 0.35-8 keV limiting unabsorbed luminosity assuming a power law with P =1.9 
and Nu =4.09 x 10^° cm-^. 



5' = 2.9 kpc 




30.0 0:55:00.0 30.0 54:00.0 

Right ascension 



Fig. 1. — An RGB rendering of our Chandra exposure of NGC 300. Soft X-ray emission (0.35-1 keV) is shown in red, medium X-ray 
emission (1-2 keV) is shown in green, and hard X-ray emission (2-8 keV) is in blue. The white outline shows the D25 isophote. 

Survey of M 33 (ChASeM33) source catalog. We present 
a summary of the approach here; the reader is referred 
to ChASeM33 work for further details. The CIAO task 
wavdetect (Freeman et al. 2002) was first used to cre- 
ate a list of potential source candidates. The positions 
of these potential sources was then used as input to 
ACIS-Extract (AE; [Broos et al. 2010), a source extrac- 
tion and characterization tool. AE was used to determine 
various source properties (source and background count 
rates, detection significances, etc.) in multiple, user- 
defined energy bands. Additionally, AE extracts spec- 
tra (with appropriate response matrices) and ex amines 
source ligh t curves for temporal variability (see Broos 



based i?-band image. For each HST field, we identified 
potential matches between the ground-based i?-band im- 
age and the F814W frames; we typically were able to 
match 3-8 stars per iJST frame. The plate solutions were 
again computed using ccmap, with r.m.s. residuals typi- 
cally less than a few hundredths of an arcsecond in both 
right ascension and declination. The Chandra field was 
searched for X-ray sources that matched the positions of 
foreground stars or background galaxies in 2MASS. Five 
sources were identified, and we obtained r.m.s. residuals 
of OVI3 in right ascension and 0V26 in declination. The 
total alignment error for each field was computed by sum- 
ming in quadrature the r.m.s residuals in the individual 
fields and the ground-based i?-band image. The results 
of our 2MASS astrometry are summarized in Table |3] 

2.3. Source Catalog Creation 

We emp loyed the same source detection strategy de- 
scribed by Tiillmann et al. (2011) in the Chandra ACTS 



et al.|2QlQ| for details). The Poisson probability ot not 
being a source {prob_no_source; hereafter pns) is provided 
by AE for each source (taking the local background un- 
certainty into account). We use these pns values as our 
source significance threshold criteria when constructing 
the NGC 300 X-ray point source catalog. 
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56:00.0 30.0 0:55:00.0 30.0 54:00.0 

Right ascension 



Fig. 2. — A POSS image of NGC 300, with the HST fields shown in red and our Chandra field of view superimposed in blue. The 
locations of our X-ray point sources are shown by green crosses. 



TABLE 3 

Alignment of Chandra and HST Images to 2MASS 



Description 


Observation ID 


R.A. 


Decl. 


# Sources Used 


r.m.s., 


r.m.s., 


# of X-ray 






(J2000) 


(J2000) 


For Alignment 


R. A." 


Decl." 


Sources 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


ground-based R 








4 


0(^00773 


0'.'0225 




Chandra 


12238 


00:54:53.05 


-37:38:42.1 


5 


0(^1332 (0(^1330) 


0(^2580 (0^2570) 


95 


Wide-1 


10915 


00:54:21.37 


-37:37:56.3 


5 


0'.'0518 (0(^0512) 


0(^0400 (0(^0327) 


5 


Wide-2 


10915 


00:54:34.7 


-37:39:25.2 


4 


0(^0171 (0(^0152) 


0(^0342 (0^0258) 


6 


Wide-3 


10915 


00:54:47.7 


-37:40:51.6 


4 


0'.'0171 (0(^0152) 


0(^0535 (0(^0485) 


5 


NGC 300-1^ 


9492 


00:55:35.66 


-37:41:25.2 


3 


0(^0340 (0(^0331) 


0(^0421 (0^0356) 


1 


NGC 300-2 


9492 


00:54:52.96 


-37:38:59.8 


6 


0'.'0772 (0(^0768) 


0(^0435 (0(^0372) 


2 


NGC 300-3 


9492 


00:54:55.33 


-37:42:04.5 


4 


0(^0179 (0(^0162) 


0(^0263 (0^0137) 


5 


NGC 300-4 


9492 


00:54:24.42 


-37:34:20.8 


5 


0(^0528 (0(^0522) 


0(^0449 (0^0388) 


2 


NGC 300-5 


9492 


00:54:26.89 


-37:39:01.4 


8 


0(^0230 (0(^0217) 


0(^0441 (0^0379) 


9 


NGC 300-6 


9492 


00:54:27.12 


-37:44:23.2 


4 


0(^0253 (0(^0241) 


0(^0400 (0^0331) 


7 



Note. — 

"Alignment errors listed show the uncertainties in each individual field and the 
error in parentheses is the uncertainty reported by ccmap for the individual field. 
^Used rxyscale. 

Our initial source list was deliberately generated with 
many more sources than we anticipated being statisti- 
cally significant. We then applied an iterative procedure 
to remove sources with increasing significance thresholds, 
until only a small number of false sources remained in 
the final source list. Each time the source list was fil- 
tered, the output was visually examined for potentially 
lost (real) sources. We used the same images described in 



ground-based i?-band image added in quadrature. The 



Section |2.1| using a variety of binning and energy ranges 
to create the initial source list: wavdetect was run on 
each energy and bin combination, using scale sizes rang- 
ing from 1 to 16 in a power-of-two sequence. To generate 
the initial list of candidate sources, we ran wavdetect 
with the following parameters: the source significance 
{sigthresh) was set to 10~^ (to catch all real sources and 
many false ones), the cleansing threshold {bkg sigthresh) 
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was set to 10 ^ ( Freeman et al.|2002 ), and we used max- 



imum number of cleansing iterations {maxiter) of 5. De- 
fault values were used for all other wavdetect param- 
eters. The resulting candidate lists were merged and 
culled of duplicate sources, resulting in an initial source 
list comprised of 799 candidate sources. 



criteria result from values below ~ 10 ^ (|Nandra et al. 


2005 


Georgakakis et al.|2008 


Tiillmannet al.|2011|). We 


adopt 


a value of pns = 4 x 10" 


as our criteria for select- 



ing those sources that are highly probable to be genuine. 
However, due to the high number of source candidates 
in our initial list, source crowding posed a serious issue 
and made defining background regions problematic. To 
address this issue, the output from our first AE iteration 
was filtered with pns = 10~^. The filtered source list 
contained revised source properties and improved back- 
ground estimates. All sources failing to meet this initial 
pns threshold were visually examined, and any obviously 
missed sources were added back to the candidate source 
list. The source list resulting from our first iteration con- 
tained 224 sources. We then ran AE again on the filtered 
source list, this time raising the pns threshold to 10~^. 
Source crowding was again reduced and spurious sources 
were filtered out, yielding an output list of 129 potential 
sources. Again, sources failing to meet our new pns cri- 
teria were visually examined and returned to the source 
list if necessary. 

AE was run a final time on the resulting source list, 
and the pns value was computed in each of the following 
nine energy bands (in keV): 0.5-8, 0.5-2, 2-8, 0.5-1, 1-2, 
2-4, 4-8, 0.35-1, and 0.35-8. To be included in the final 
NGC 300 catalog, a sources was required to have a pns 
value less than 4 x 10~^ in any of the nine energy bands; 
if only the 0.35-8 keV band were considered, ~4% of sig- 
nificant sources would have been lost. The final CLVS 
source catalog for NGC 300 contains 95 source^ Their 
positions and properties, such as detection significance 
{pns value), net counts, and photon flux (in eight differ- 
ent energy bands) are listed in Tables pjlTl 



2.4. Sensitivity Maps 

AE provides an estimate of the photon flux {flux2) based 
on the net number of source counts {net_cnts)^ the expo- 
sure time, and the mean ARF in the given energy band. 
To convert this photon flux into an energy flux, knowl- 
edge of the intrinsic X-ray spectrum is required - how- 
ever, less than 20% of the NGC 300 X-ray sources have 
a sufficiently large number of counts that can be used to 
constrain a spectral model. 

To enable a uniform comparison between NGC 300 
and M33, we ad o pt th e same approach as employed by 
Tiillmann et al. (2011): we assume an absorbed power 



law spectrum, with F = 1.9 and = 4.09x10^^ cm~^, 
to estimate the energy flux for each source. This model 
is hereafter referred to as the "standard model," and ap- 
pears to be reasonable for both XRBs associated with 
NGC 300 as well as for background AGN (especially at 
low flux levels). We note that the energy flux estimate 
for individual sources may be off by as much as a factor 
of two if the true spectrum is different from our assumed 

^ Source catalogs are available in FITS format at: ,http : //www. | 
[astro . W ashington . edu/ users/bbinder/ CLVS/ 



one, although the net effect of this assumption is not ex- 
pected to systematically bias the logTV-logiS relation to 
higher or lower fluxes. 

The following analyses require, for each point in our 
Chandra exposure, a sensitivity map which provides the 
energy flux level at which a source would be detectable. 
The CIAO task liin_sens was used to compute limiting 
sensitivity maps from the background image, exposure 
map, and ACIS-I PSF image. Our 90% complete lim- 
iting luminosities in the 0.35-8 keV, 0.5-2 keV, and 2-8 



keV ranges are 1.5x10"^^ erg s~-^, 8.8x10"^^ erg s~^, and 
4.1x10^^ erg s~^, respectively. Table N summarizes the 
limiting luminosities for the 70%, 90%, and 95% com- 
pleteness levels for our total energy range, the soft band, 
and the hard band. 

3. THE X-RAY SOURCE CATALOG 

Various properties of individual X-ray sources, such 
as their temporal variability or spectral shape, can be 
used to constrain the physical origin of the X-ray emis- 
sion. Additionally, bulk population properties (i.e., ra- 
dial source distribution and the XLF) allow us to discrim- 
inate between those sources associated with NGC 300 
and background AGN in a statistical sense. We expect 
the observable X-ray point source population of NGC 300 
to be comprised of HMXBs, LMXBs, background AGN, 
and a few bright supernova remnants (SNRs). To match 
each point source to its most likely physical origin, we 
use our catalog to identify time- variable sources, analyze 
hardness ratios (for faint sources) and spectra (for suffi- 
ciently bright sources), and create XLFs. The additional 
matched HST optical observations allow us to identify 
candidate optical counterparts and measure correlations 
between the X-ray emission and the underlying stellar 
age as a function of radius of NGC 300. 

3.1. Comparison with XMM-Newton and Temporal 
Variability 

Some classes of XRBs, as well as background AGN, ex- 
hibit long-term X-ray variability on timescales of months 
to years, while SNRs and some HMXBs will show persis- 
tent X-ray emission. To search for long-term variability, 
we compared our Chandra X-ray catalog to the XMM- 
Newton survey conducted by ^Carpano et al.| ( |20Q5| , 
which contained 86 sources within the 1^25 isophote of 
NGC 300. We consider sources with X-ray positional 
offsets < 3'.' 5 to be matches. 

To correct for the differences in energy ranges between 
Chandra (0.35-8 keV) and XMM-Newton (0.3-6 keV), we 
converted all XMM-Newton fluxes into Chandra fluxes 
assuming our standard model. Due to the differences 
in energy resolution, detector efficiencies, etc., we can- 
not compare Chandra and XMM-Newton hardness ra- 
tios to search for changes in spectral shape (see Sec- 
tion |3^ fo£further_discussion of source hardness ratios). 



Tiillmann et al. (2011) defined a variability threshold as 
77 = (tluxmax - tluXminj/Aflux, whcrc the flux error Afiux 
is calculated using the Gehrels approximation (appropri- 
ate for low count data; Gehrels 1986), to identify sources 
with potential long-term variability. Any source with 
7^ > 5 is considered potentially variable; according to 
this criteria, 8 sources in the NGC 300 catalog exhibit 
long-term variability. Table [9] lists the temporal variabil- 
ity measurements for the 52 sources observed by both 
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TABLE 4 
NGC 300 Source List 



Source 


Source ID 


R.A. (J2000) 


Decl. (J2000) 


Positional 


Total exp. map 


RsTC 





No. 




n 


n 


Error (") 


value (s cm^) 


(sky pixel) 


n 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


1 


005411.98-373951.8 


13.549937 


-37.664408 


0.70 


1.589x10'^ 


13.43 


8.2 


2 


005412.31-373359.6 


13.551319 


-37.566583 


0.93 


1.434x10'^ 


16.66 


9.3 


3 


005413.98-373710.5 


13.558268 


-37.619588 


0.36 


1.588x10'^ 


12.14 


7.9 


4 


005415.57-373316.4 


13.564907 


-37.554563 


0.59 


1.210x10'^ 


15.25 


9.2 


5 


005419.92-373744.5 


13.583001 


-37.629032 


0.61 


1.711x10'^ 


8.49 


6.6 



Note. — Column 2: the source ID also contains the source coordinates (J2000.0). Column 5: the positional uncertainty is a simple 
error circle (in around R.A. and decl. Column 6: sum of mean exposure map values in the source region. Column 7: the average radius 
of the source extraction region (1 sky pixel = 0(^392). Column 8: off-axis angle. 
Only the first five entries are shown. 



TABLE 5 

Logarithmic pns Values for Different Energy Bands 



Source 
No. 
(1) 


\og(pns [1]) 
(0.5-8.0 keV) 
(2) 


\og(pns [2]) 
(0.5-2.0 keV) 
(3) 


\og(pns [3]) 
(2.0-8.0 keV) 
(4) 


\og(pns [4]) 
(0.35-8.0 keV) 
(5) 


\og(pns [5]) 
(0.35-1.1 keV) 
(6) 


\og(pns [6]) 
(1.1-2.6 keV) 
(7) 


\og(pns [7]) 
(2.6-8.0 keV) 
(8) 


\og(pns [8]) 
(0.35-2.0 keV) 
(9) 


1 


-9.32 


-8.60 


-3.16 


-9.13 


<-10 


<-10 


-2.05 


-2.05 


2 


-4.54 


-3.57 


-2.17 


-4.81 


-0.42 


-2.80 


-2.77 


-2.77 


3 


<-10 


<-10 


<-10 


<-10 


<-10 


<-10 


-6.03 


-6.03 


4 


<-10 


<-10 


-6.94 


<-10 


-3.46 


<-10 


-4.76 


-4.76 


5 


-5.48 


-4.58 


-2.15 


-6.04 


-3.20 


-3.49 


-1.41 


-1.41 



Note. — The (logarithmic) pns value is the Poisson probability of not being a source, hog(pns) values smaller than -10 have been 
replaced by <-10. All these cases are highly significant detections. 
Only the first five entries are shown. 



TABLE 6 

Total Net Counts in Different Energy Bands 



Source net-cnts[l] 
No. (0.5-8.0 keV) 
(1) (2) 



22.9^ 

7 q+6.5 

•'^-5.4 

:,+8.5 
-7.4 

11 2+^-^ 



net-cnts[2] 
(0.5-2.0 keV) 
^ 



13.3-^ 

o 0+4.5 

45.6l^-? 
25.5 



net-cnts[3] 
(2.0-8.0 keV) 
W 



+5.U 
-3.9 
+ 5.4 
4.3 
+ 6.6 
5.5 
^+6.3 
^-5.2 
+ 3.8 
2.6 



net-cnts[4] 
(0.35-8.0 keV) 



+6.6 
5.5 
+6.6 
5.5 
+9.8 
-8.8 
+8.5 
7.4 
+5.0 
3.8 



net-cnts[5] 
(0.35-1.1 keV) 
(6) 



net-cnts[6] 
(1.1-2.6 keV) 

m_ 



1 7 1 +^-^ 
-^'•-^-4.3 

-+6.2 
-5.1 



net-cnts[7] 
(2.6-8.0 keV) 
W 



/1+4.5 
-3.3 
□+5.3 
-4.1 

14.011^ 

:; + 5.7 
-4.6 
9+3.4 
-2.2 



net-cnts[8] 
(0.35-2.0 keV) 
W 



13.2-^ 
45.6i^:^3 

,+6.3 
-5.2 



70.9: 
45.3^ 



6.3 



+6.3 
5.2 
+3.8 
2.6 



9 

9.1 
25.3 
19 
4.9; 



22 

18.2 

70 

45.1 

12.1 



< 1.1 

+2.7 
1.3 
+5.2 
4.1 
+3.8 
2.6 
+3.2 
1.9 



0.7; 

16.5 
6.0 

3.7; 



24.6^ 



5.2 



+3.6 
-2.4 



6.4: 

9.9^ 

14.0; 

14.5: 
3.2j 



25.3^ 



7.2- 



+4.0 
-2.8 



Note. — Only the first five entries are shown. 



TABLE 7 

Photon Flux (photons cm~^ s) in Different Energy Bands 



Source 


\og{flux[l]) 


\og{flux[2]) 


\og{flux[3]) 


\og{flux[4]) 


\og{flux[b]) 


\og{flux[6]) 


\og{flux[7]) 


\og{flux[8]) 


No. 


(0.5-8.0 keV) 


(0.5-2.0 keV) 


(2.0-8.0 keV) 


(0.35-8.0 keV) 


(0.35-1.1 keV) 


(1.1-2.6 keV) 


(2.6-8.0 keV) 


(0.35-2.0 keV) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


1 


-5.75 


-6.11 


-6.09 


-5.74 




-6.07 


-6.25 


-6.07 


2 


-5.84 


-6.28 


-6.08 


-5.81 


-6.83 


-6.40 


-6.04 


-6.20 


3 


-5.26 


-5.57 


-5.67 


-5.25 


-5.54 


-5.69 


-5.91 


-5.54 


4 


-5.28 


-5.65 


-5.60 


-5.27 


-5.77 


-5.74 


-5.73 


-5.61 


5 


-6.07 


-6.45 


-6.39 


-6.03 


-6.23 


-6.60 


-6.57 


-6.35 



Note. — Only the first five entries are shown. 

Chandra and XMM-Newton^ with sources showing evi- 
dence of variabiUty hsted in boldface. 

We found a total of 75 X-ray sources that were de- 
tected by either Chandra or XMM-Newton, but not both. 
For these sources, we use the detection limit of the ob- 
servatory in which the source was not observed to set 
a lower-limit on the long-term source variability. X-ray 
transients (XRTs) observed in the Milky Way and other 



nearby galaxies originate in either LMXBs or HMXBs, 
but in all cases the peak X-ray luminosity exceeds the 
quiescent X-ray luminosity by more than an order of 
magnitude. We therefore define XRT candidates in 
NGC 300 as those sources whose maximum luminosity is 
at least an order of magnitude above the detection limit 
for the exposure in which they were not observed. The 



25 sources meeting this criteria are listed in Table 10 
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TABLE 8 

Limiting Luminosities for our NGC 300 Observation 



Energy 


70% 


90% 


95% 


(keV) 


(1035 erg s-i) 


(10^5 erg s-i) 


(10^5 erg s-i) 


(1) 


(2) 


(3) 


(4) 


0.35-8 


7.7 


15 


45 


0.5-2 


4.5 


8.8 


27 


2-8 


21 


41 


120 



TABLE 9 

NGC 300 Cross Identifications with XMM-Newton Sources 



CLVS 


XMM-Newton 


Positional 




T] 


fluXmax/fluXmin 


Source No. 


Source No. ID 


Oflfset {") 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


1 


41 


0.76 


0.7958 


0.21 


1.41 


3 


17 


1.38 


0.9152 


0.79 


2.49 


6 


45 


0.98 


0.5160 


0.44 


1.69 


7 


39 


0.60 


0.5590 


2.79 


4.30 


8 


42 


1.70 


0.2267 


0.14 


1.26 



Note. — Columns 1 and 2 give the CLVS and XMM-Newton source numbers, respectively. Column 3: the difference (in arcseconds) 
between the CLVS and the XMM-Newton source positions. Sources with positional offsets less than 3^5 are considered matches. Column 4 
is the K-S probability of the source being constant within our Chandra observation (short-term variability). Sources with K-S probabilities 
of £ < 3.7 X 10~^ are likelv to show short-term variability and are shown in boldface. Column 5: 77 is the variability index defined by 
[Tiillmann et al.| ( |2011| > and is sensitive to long-term variability; sources with variability indices of 77 > 5 are likely to show long-term 
variability and are sliown in boldface. Column 6: variability factor expressed as the ratio of fluxmax/fluxmin- 
Only the first five entries are shown. 



We additionally searched the Chandra source catalog 
for short-term variability in the 0.35-8 keV band (i.e., 
variability that occurs on timescales less than the expo- 
sure time, ~18 hours). We used a Kolmogorov-Smirnov 
(K-S) test to compare the cumulative photon arrival time 
distribution for each X-ray source to a uniform count rate 
model. The K-S test returns a probability that both 
distributions were drawn from the same underlying dis- 
tribution. One one source (18) has a K-S probability 
indicating rapid variability at the 3a level. 

3.2. Hardness Ratios 

The X-ray spectrum of a source can provide a key di- 
agnostic for separating different sources. For example, 
XRBs and AGN are typically well described with power 
laws with photon indices T =1.5-2.0, while soft X-ray 
sources (such as SNRs) show spectra consistent with a 
thernial pl asma of temperatures of ~ 0.3 — 0.8 keV (Long 
et al. |2010p . Still other sources, such as low luminosity 
AGN , may exhibit an X -ray spe ctrum consistent with 
bremsstrahlung emission (Yi & Boughn 1998). However 



in order to constrain the fit parameters of these various 
models to believable values, a sufficiently high signal- 
to-noise is required. Even in the best case scenario, low 
signal-to-noise data (less than ~50 counts) can only have 
a power law photon index constrained to ~50% accuracy. 
The majority of the NGC 300 X-ray sources in our cat- 
alog fall into this low-count category (see Figure [3|. 

To assess the X-ray spectral shape for low-count 
sources, a different approach must be used. Hardness 
ratios (HRs), also referred to as X-ray colors, are com- 
puted for various energy bands to measure the fraction 
of soft {S, 0.35-1.1 keV), medium (M, 1.1-2.6 keV), and 
hard (iJ, 2.6-8.0 keV) photons emitted by a source. A 
source comprised of a thermally emitted X-rays, such 
as a SNR, will show characteristically soft X-ray colors, 
while an obscured background AGN (having had most 



of its softer X-rays absorbed) will exhibit a significantly 
harder X-ray color. 

We evaluate two HRs for each source using source and 
background counts determined by AE in the bands listed 
above. We define a "soft" color. 



HRl 



and a "hard" color. 



HR2 



M-S 

WTWTs' 



H-M 

WTWTs' 



(1) 



(2) 



We use 
(BEHR, 



the 



Bayesia n Estimation of Hardness Ratios 
in an approach described by 

nlrTse< 



Park et al. 



Tiillmann et al.| ( |2Ull ) "(see also iPrestwich et al.||2009|) 

to determine the hardness ratios of the NGC 300 X-ray 



ray 

sources. The BEHR code has the advantage of requir- 
ing source and background counts to be non-negative 
(net counts for low-count sources can be negative due 
to fluctuations in the source and background estimates, 
resulting in unphysical HRs); however, it cannot directly 
compute HRs in the form we require. We therefore use 
the BEHR code to generate 50,000 samples of the proba- 
bility distribution of source fluxes for each energy band. 
To do this, we use as input to BEHR the following quanti- 
ties from the AE output for each source: source and back- 
ground counts, the AE "backscale" parameter (to account 
for the ratio of the source and background extraction ar- 
eas and efficiencies), and a factor converting from counts 
to flux (i.e., the exposure time multiplied by the mean 
of the ARE over the extraction region in the appropriate 
energy band). We then run the BEHR code twice, once 
for S and M and again for S and H. 

We combine the results from our two BEHR runs to 
obtain 50,000 values of HRl and HR2 for each source. 
Because the difference of two Poisson distributions does 
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TABLE 10 

Candidate Transient X-ray Sources from Chandra and XMM-Newton 





Xv.-Ti.. I JZUUU^ 


LJkH^L. ^jzuuu^ 


T] 


fliiY /fliiY ■ ^ 

nuXmax/ nuXnim 








(A^ 


\^) 






-0 ( .OO^ODO 


1 1 1 

ID. 1 1 


91 7fi 
Zl . ( D 




iO.OUOo ( ± 


-0 ( .04t04tUU 


I 1 1^9 

I I .OZ 


1 ^ 9.A 




±0.04t±oOO 


-0 ( . ( yoooo 


±0.0 ( 


9/L 5^1 


r^vo 9Q 


io.oo ( lyo 


07 K/lOOQ/l 

-0 ( .o^zzy^ 


0/1 

y .y4t 


1 9 i=iO 
IZ.OU 




lO.DOyoiO 


Q7 c;Q/i/if^o 

-0 ( .OO^^KJO 


^^9 9*^ 


/19 i^;*^ 
^z. 00 




io.ooyozu 


Q7 c:/197QC=; 

-0 ( .o4tz ( 00 


19 01 
iz.y 1 


1 S ^A 


pxn f^9 


1 7074^=17 


-0 ( .0100^ ( 


12 24 


1 7 9fi 

1 ( .ZD 


V j\ V y - ( o 


1 770901 
10. ( ( UZUl 


-0 ( .oyo'iD^ 


1 40 
10. ^y 


90 /17 


V y\ V y - ( D 


1 77009"=; 


Q7 c^DzlOl 7 

-0 ( .ou^y 1 ( 


0.00 


1 fi*^ 

lU.DO 


pvo 07 
LyyvW-o ( 


10.oZ04t4tO 


Q7 c;c;qo79 
-0 ( .oooy ( z 


1 9 09 

iz.yz 


1 S 5^0 
lo.oU 




10.00 ( ouo 


-0 ( .^0000^ 


41 25 


Do. 00 


yvlVllVl-Z 


1 70'=;9t=;o 
10. ( yozou 


-0 ( .ouyooD 


10. OU 


^UO.OO 


yvlVllVl-O 


1 70070S 
10. ( uy ( uo 


-0 ( .D^ ( UoO 


91 OR 
Zl .UO 


oyi.oy 


vix/rix/r A 

^lVllVl-4t 


±0. ( yo4too 


Q7 «/La/L1 7 


11 11^ 
11.10 


9.'^ 9f^ 
OO.ZO 




1 fifiOzLI 7 
lo.uuy^i ( 


Q7 fi»09f=;0 
-0 ( .UoUZOU 


^ "^0 
o.oy 


f=;7 7fi 


XMM-8 


13.835292 


-37.802833 


3.86 


35.67 


XMM-10 


13.835292 


-37.643944 


6.35 


60.56 


XMM-16 


13.706458 


-37.666861 


5.55 


51.68 


XMM-18 


13.700292 


-37.782611 


5.00 


43.24 


XMM-23 


13.583208 


-37.652472 


4.15 


30.07 


XMM-25 


13.926000 


-37.673250 


5.91 


28.49 


XMM-43 


13.722000 


-37.722000 


3.27 


18.90 


XMM-52 


13.566292 


-37.566292 


2.92 


12.68 


XMM-53 


13.807708 


-37.807708 


3.27 


12.51 


XMM-59 


13.849958 


-37.756694 


3.28 


12.29 



Note. — - ^Sources labeled as 'CXO' were detected in our 
|et al.| ( |2005t XMM-Newton catalog only. 

^The long-term variability index 77 and the variability factor 



Chandra catalog only; sources labeled as 'XMM' were detected in the |Carpano| 
fluxmax/fluxmin Hstcd are lower limits. 




0.50 
0.40 
0.30 
0.20 
0.10 

000 . ....... 

50 100 150 

Net Source Counts (0.35-8 keV) 

Fig. 3. — Histogram of the net counts observed in the 0.35-8 keV band for our NGC 300 X-ray catalog, excluding source 80 (NGC 300 
X-1). We find the majority of our sources have less than 50 counts, making spectral modeling for these sources unreliable. 

not follow Poisson statistics, we do not use the error es- 
timates directly provided by these BEHR runs. Instead, 
the HR value is defined as the mean of the distribution, 
and the credible interval is evaluated based on the 68.2% 
equal-tail estimates (i.e., 0.682/2 of the samples have val- 
ues below the lower limit, and 0.682/2 of the samples 
have values above the upper limit). Our HRs for each 
source are listed in Table [TTl 

To aid in the interpretation of our HR calculations, we 
employ an X-ray color-color classification scheme simi- 
lar to that used in [Kilgard et al.| ([2005|. Due to the 



differences in energy bands considered in this work, we 
modify the classification scheme slightly. We simulated 
spectra for thermal sources of varying temperature (0.1-2 
keV) and power law sources (with T =0.4-3) with differ- 
ent levels of absorption and calculated HRl and HR2 for 
each source to inform our X-ray color cuts. We define 
six categories: XRBs (which likely also contain a sig- 
nificant fraction of background AGN), SNRs, absorbed 
sources ('ABS'), indeterminate soft sources ('SOFT'), in- 
determinate hard sources ('HARD'), and sources with an 
indeterminate spectral type (labeled 'INDET', although 
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TABLE 11 
X-RAY Source Hardness Ratios 



Source No. 


HRl 


HR2 


Class 


(1) 


(2) 


(3) 


(4) 


1 


0.72 ± 0.25 


-0.71 ± 0.32 


ABS 


2 


0.19 ± 0.41 


-0.09 ± 0.51 


XRB 


3 


0.36 ± 0.13 


-0.56 ± 0.12 


SOFT 


4 


0.50 ± 0.15 


-0.47 ± 0.21 


ABS 


5 


-0.07 ± 0.46 


-0.15 ± 0.35 


XRB 



Note. — Only the first five entries are shown. 

we do not observe any sources with hardness ratios fitting 
this criteria). Indeterminate sources are those whose X- 
ray spectra are ambiguous; for example, 'HARD' sources 
may have exceptionally flat power law spectra, suffer 
extreme absorption, have multiple spectral components 
that cannot be distinguished due to our low number of 
source counts, or some combination of these effects. The 
color-color classification scheme is summarized in Ta- 
ble [I2l and a plot of the HRs for the NGC 300 X-ray 
source population is shown in Figure |4] 

We find most of the NGC 300 sources fall within 
the 'XRB' category (41%), and an additional ^12% of 
sources are categorized as SNRs. This is not surprising - 
X-ray emission from non-active spiral galaxy with recent 
star formation is expected to be dominated by XRBs and 
bright SNRs. Nearly one-third of the NGC 300 X-ray 
sources show evidence for absorption beyond the Galac- 
tic column. Two sources exhibit unusually hard X-ray 
colors, and are likely background AGN whose soft X-ray 
emission has suffered a high degree of absorption. Thir- 
teen sources show unusually soft X-ray colors. These may 
be very soft XRBs, sources with unusual absorption, or 
foreground objects not associated with NGC 300. 

3.3. X-ray Spectral Analysis 

As discussed above, the X-ray spectrum of a source 
depends on the physical X-ray production mechanism. 
The value of the spectral fit parameters can reveal further 
clues to the nature of the X-ray emission. For example, 
XRBs containing a NS primary will show a hard power 
law with r ~ 1 — 1.5, while BH primaries produce a 
significantly softer spectrum with F ~ 2 — 2.5. 

AE utilizes the CIAO tools dmextract, mkacisrmf , and 
mkarf when extracting spectra and generating response 
products. Due to the low number of counts found for 
most of the sources in our catalog, we fit spectral models 
to the unbinned spectra and use C-statistics in lieu of 
traditional statistics. The only X-ray source to have 
more than 200 net counts in the 0.35-8 keV band is source 
80, the HMXB NGC 300 X-1. Detailed spectral fitting 
and timing analysis was carried out for this source in 
[Binder et al.'('2011aV 

The ungrouped spectra for 53 of the 95 sources were 
automatically fit in AE with one of four models : a power- 
law model, an APEC model ( Smith et al.||200l' ), or either 
a power-law or thermal model with additi onal absorp- 
tion. The tbabs model (Wil ms et al.p OOO*) was used to 
model the Galactic absorption (A^h =4.09x10^^ cm~^; 



Kalberla et al. 2005). To determine the best-fit model, 
we first rejected those models which predicted unreason- 
ably high photon indices (F > 4) or plasma temperatures 
{kT > 6 keV). To assign a goodness-of-fit comparison for 
each of the models using C-statistics, we used XSPEC's 



goodness command to perform 5,000 Monte Carlo cal- 
culations of the goodness-of-fit for each model. The com- 
mand simulated spectra based on the model and com- 
puted the fraction of simulations with a fit statistic less 
than that found for the data - if the observed spectrum 
was produced by that model, the goodness should be 
around 50%. We therefore chose the best-fit model as the 
one for which the goodness command returned closest 
to 50%. In some cases, the goodness command returned 
values near 50% for multiple models, or did not return 
a reasonable value for any model. In these cases, we 
checked each spectrum for indications for emission fea- 
tures. In cases where evidence of thermal emission was 
present, the thermal model (with fixed A^h) was adopted; 
for all other sources, the power law with fixed A'h was 
chosen. The best-fit model and parameters for sources 
in NGC 300 with greater than 50 net counts are listed in 
Table [m 

3.4. Radial Source Distribution 

We assign an inclination-corrected galactocentric dis- 
tance to each X-ray sources assuming th e following 
galaxy parameters for NGC 300, as in Gog arten et al. 
|2QTq[ ): ao = 13?722833, So = -37?684389 (galaxy cen- 
ter), i = 42°, (in clination), and = 111° (position angle 
Kim et al. 2004| ). The maximum radius we find for an 
X-ray source m our Chandra field is r = 8.06 kpc. 

We divide the X-ray sources into radial bins based on 
their inclination-corrected distance from the center of 
the galaxy. The bins are spaced 0.9 kpc apart in ra- 
dius so that we can directly compare the X-ray radial 
di strib ution to the stellar population analysis presented 
in Gogarten et al. (201 0J. We use the cumu lative logA^- 
logS distribution of Cappelluti et al. (2009) to estimate 
the expected contamination by background AGN, which 
should be nearly flat with galactocentric radius. 

The radial X-ray source distribution is shown in Fig- 
ure |5j and clearly shows that the number of X-ray sources 
per square kiloparsec decreases with galactocentric ra- 
dius. Since NGC 300 does not possess a substantial 
bulge, we fit the background-subtracted radial X-ray 
source distribution to a simple exponential profile, and 
find the best-fit relation to be 



^x{r) = (12.8 ± 3.0)e-"/(^-^±°-') - (0.57 ± 0.85), (3) 

where A/x (r) is the number of X-ray sources per kpc^ ex- 
pected at a given galactocentric radius, r (in kpc). The 
exponential scale length is 1.7± 0.2 kpc, similar to the 



1.3±0.1 kpc scale length found by Gogarten et al. (2010) 
and in excellent agreement with estimates derivea trorn 
the i^T-band surface brightness ( ^Muhoz-Mateos et al 
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TABLE 12 

Hardness Ratio Source Classification Scheme 



Classification 


Definition 


# Sources 


(1) 


(2) 


(3) 


X-ray binary ('XRB') 


-0.4 < HR2 < 0.4, -0.4 < HRl < 0.4 


39 


Absorbed source ('ABS') 


HRl > 0.4 


30 


Supernova remnant ('SNR') 


HR2 < 0.4, HRl < -0.4 


11 


Indeterminate hard source ('HARD') 


HR2 > 0.4, -0.4 < HRl < 0.4 


2 


Indeterminate soft source ('SOFT') 


HR2 < -0.4, -0.4 < HRl < 0.4 


13 



1.0 



0.5 



^ 0.0 

II 



X -0.5 



-1.0 

-1.0 -0.5 0.0 0.5 1.0 

HR2 = (H-M)/(S+M+H) 




Fig. 4. — X-ray color-color diagrams of all X-ray sources detected in NGC 300. The thick cross in the upper-right corner shows the 
typical size of the errors. The diagram has been broken into preliminary source identification regions, described in Table [12] The regions 
are labeled as: 'ABS' for absorbed sources, 'XRB' for X-ray binaries, 'SNR' for supernova remnants, 'SOFT' for indeterminate soft sources, 
'HARD' for indeterminate hard sources, and 'INDET' for sources with an indeterminate spectral shape (although we do not observe any 
sources in this category). The red lines indicate the zero-count limits in the hard ('H'), medium ('M'), and soft ('S') bands. 

ipated to dominate the X-ray sources at greater than ~4 
kpc, the 'XRB' class hkely suffers from AGN contamina- 
tion. 'SOFT' sources additionahy span a wide range in 
radii, with a median distance of 4 kpc, suggesting that 
at least some of these sources are likely to be associ- 
ated with NGC 300. Absorbed and indeterminate hard 
sources have median galactocentric radii greater than 4 
kpc, implying that many of these sources are background 
or foreground objects. 



2007). We additionally find the X-ray scale length to 
be consistent with multiwavelength measurements. The 
radial source distribution suggests that the X-ray sources 
in NGC 300 are the evolutionary end-points of a parent 
population of disk stars. 

For an exponential disk, we expect 90% of stars to 
be inside of 4 kpc (6^.9), suggesting that the majority of 
sources beyond this radius are not associated with the 
stehar disk of NGC 300. This suggests that 37 X-ray 
sources beyond 4 kpc are likely to be AGN (i.e., 40% 
of the sample). This interpretation is consistent with 
the estimated number of background and hard sources 
identified in our hardness ratio analysis. 

We use the source classifications determined by our 
hardness ratio analysis to investigate the differences in 
galactocentric radii among the different types of sources. 
The results are summarized in Table [l4j We find both 
the preliminary SNR and XRB candidates to have me- 
dian radii less than 4 kpc, suggesting that many of these 
are indeed associated with NGC 300. However, XRB- 
classified sources span a wide range of galactocentric radii 
(out to 6.7 kpc). Given that background AGN are antic- 



4. SOURCE CLASSIFICATION 
4.1. Supernova Remnants 



In Section |3.2[ we found that 11 X-ray sources in 
our catalog had hardness ratios consistent with those of 
SNRs. Of these, only one source (#43) had a sufficient 
number of counts for spectral modeling to be performed. 
Source 43 was best-fit by a thermal plasma model with a 
tempera ture kT=OA keV, consistent with SNRs observed 
in M33 ( Long et al.||2010 ). Taking this as our "standard 
SNR" model, we estimate the unabsorbed 0.5-2 keV lu- 
minosities for the remaining 10 SNR candidates and find 
a range in luminosities of ~ 5 x 10^^ lum to ~ 6 x 10^^ 
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TABLE 13 

Best-Fit Spectral Models for Sources With More Than 50 Counts 



No. 


Best-fit^ 
ivioaei 


(cm j 


r 


(^Ke V ) 


Norm^ 




dof^ 


C 


Lx (0.5-2.0 keV)^ 
(A n36 


Lx (0.5-8.0 keV)-^' 
erg s~^) 


"goodness^" 
[/o) 
















\ ' ) 










o 
o 


pow 




1 7-1-0 




O. IZ X lU 


6 


±u±o 


^7f\ 1 Q 


o.ozCi. ± 


^-+3.3 


40.00 




pow 




1 f^-l-fl A 




9 QQ V 1 n~ 


6 


lUlU 


UZiJ.O ( 


9 «±n 7 

Z.U U. ( 


Q+2.2 


oZ.oU 


23 


pow 


n 7+0.6 


9 9+0.9 




l.lOx 10~ 


5 


1016 


464.43 


9 7+1.1 
^- ' -0.9 


10.3_3 3 


32.70 


24 


£ipGC 






1.22x10" 


5 


1016 


750.98 


3 4+2-? 


7 9+i-f 

' -^-3.8 


98.38 


25 


pow 


<0.4 




4.42x 10~ 


6 


1016 




q 9+1.9 


n 9 + 3.8 


55.28 


Zo 


pow 




1 A -Un o 

i.4ztU.Z 




Z. i < X iU 


6 


iUio 


/I QQ 

4oo.DU 


9 /I -Un 7 


9+2.7 


A A KA 
44.04 


33 


pow 


2+2-^ 


1 8+2-? 




6.99x10" 


6 


1015 


679.83 




16 7+?;?. 
' -10.5 


55.04 


35 


pow 


1 1+^-° 

— 0.5 


l-3+oI 

— U.4 




6.66x10- 


6 


1015 


706.75 


7.4+1-0 

— z.o 


27. 5 1 n r> 


41.04 


39 


pow 




1.5±0.4 




2.29x10- 


6 


1016 


590.38 


^•^-0.5 


7 7+2.5 
' • ' -1.8 


64.26 


43 


apec 






0.4±0.1 


1.92x10- 


5 


1016 


456.08 


6.8+1-27 




66.84 


50 


pow 




1.6±0.3 




4.12x10- 


6 


1016 


631.44 


4 4+0-8 


12.3+1? 


91.80 


52 


pow 




2.51M 




8.22x10- 


6 


1016 


390.55 




5.711.1 


55.88 


53 


pow 


<0.3 






3.31x10- 


6 


1016 


380.30 


^+1.7 




52.93 


87 


apec 


< 0.2 




2 6+^-^ 

^•"-0.8 


7.31x10- 


6 


1016 


371.56 


9 7+1.6 
^- ' -0.8 




72.86 


89 


pow 




'-'•'^-0.3 




1.62x10- 


5 


1015 


564.66 


' "^-12.0 


97 7+7.0 
^'•'-12.7 


71.14 


95 


apec 








1.79x10- 


6 


1016 


541.01 


7+0.3 
^- ' -0.6 


7+0.5 
^- ' -0.6 


98.80 



Note. — 

The best-fit model is either a powerlaw or vapec model, depending on which one has the lower C/dof. 
Intrinsic source absorption (in units of 10^^ cm~^), if beyond the Galactic column was required. 



^ Normalization constant of the fit, in units of photons keV cm 
respectively. 

^ Degrees of freedom. 

^'•^ Unabsorbed X-ray luminosities, assuming a distance to NGC 300 of 2.0 Mpc ( Dalcanton et al.|2009 |. 
^ Results of the XSPEC "goodness" command, run using 5000 realizations. 



and cm ^ for the power law and thermal plasma models. 



TABLE 14 

Galactocentric Radii By Hardness Ratio IDs 



Source classification 
(1) 


Tmin (kpc) 

(2) 


rmax (kpc) 

(3) 


rav (kpc) 
(4) 


SNR 


0.4 


6.0 


2.5 


XRB 


1.3 


6.7 


3.4 


ABS 


1.2 


8.1 


4.1 


SOFT 


2.1 


6.2 


4.0 


HARD 


2.8 


6.9 


4.9 



erg s~^, also consistent with observations of M33 SNRs. 

To see if any of our SNR candidates have been previ- 
ously identified, we cross-correlated our 11 sources with 
the n i ultiwa vele ngth NGC 300 SNR catalogs of Bl air fc] 
Long (1997) and Pannuti et al. ([2QQQ|) an d the catalog of 
NUU3imT II regions of |Deharveng et al [ ( [r988| ). We find 
7 of our 11 sources are located within 3.5 of a SNR or 
H II region, and two of our X-ray sources are coincident 
with optically confirmed SNRs. Table 15 summarizes our 
SNR candidates. Figure [g] shows 0^5x0.5 postage stamp 
images of the Ha emission coincident with each X-ray 
source, with the soft 0.35-1 keV X-ray contours superim- 
posed. 

4.2. Optical Counterparts 

We searched for optical counterpart candidates for each 
X-ray source that was contained within one or more 
overlapping HST field. For many of the X-ray sources, 
we found multiple optical counterpart candidates within 
the Chandra X-ray error circle (typically on the order of 
~0'.'5); in these cases, the optical counterpart candidates 
are labeled 'a,' 'b,' 'c,' etc. in order of increasing dis- 



tance away from the X-ray source location. We kept only 
those optical sources that met the q uality cu ts described 
in (William s et al.| (|2009) and Gogarten et al. (2010) (i.e., 
sources with sutticient signal-to- noise, not flagged as un- 
usable, meeting predetermined crowding and sharpness 
thresholds, etc.). 

For each source, we compute a logarithmic X-ray-to- 
optical fiux ratio log{fx / fv)i where fx is the 2-10 keV 
X-ray fiux and and the fy is the optical fiux measured 
in either F555iy or F606iy, whichever is available. Al- 
though the F606VK has a slightly larger bandpass than 
the FbbbW filter, we estimate that for the majority 
of stars accessible in our HST exposures the affect on 
logjfx / fv) will be small, ~ 0.1. 

Table [l6]lists the optical and X-ray properties of all the 
X-ray source candidate optical counterparts. Figure [7| 
shows RGB-rendered b" x finding charts for each 01 
the 32 X-ray sources with overlapping HST fields. The 
X-ray source position is shown by a white cross, with a 
white circle indicating the extent of the Chandra error 
circle. Red lines are used to indicate the likely optical 
counterpart; in cases of multiple optical counterparts, the 
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Fig. 5. — The galactocentric source density profile of NGC 300 in the 0.35-8.0 keV band. The blue histogram shows the total number 
of sources per square kiloparsec. The red dashed line shows the anticipated number of background AGNs, while the black histogram 
shows the expected number of sources after correcting for background AGN. The thick dot-dashed line shows our exponential fit to the 
AGN-subtracted radial source distribution, with the lighter dot-dashed lines indicating the la uncertainty. 



TABLE 15 

Candidate Supernova Remnants and H II Regions 



Source No. 


Counterpart 


Separation 


0.5-2 keV unabs. 


Notes 








luminosity (erg s"-*^) 




(1) 


(2) 


(3) 


(4) 


(5) 



7 

15 
32 

34 

41 
43 
46 
54 
77 
81 

85 



SNR 
HII 

SNR 
HII 

SNR 

SNR 

HII 
SNR 
H II 



^6 
3'/2 

3'/3 

c/e 

2'/8 

3'/0 

3'/l 
0'.'3 

3(^4 



6.9xl035 
6.4x10^5 
1.1x10^^ 
1.5x10^6 

1.5x10^6 
5.6x10^^ 
6.3x10^^ 
1.8x10^6 
7.5x10^5 
4.7x10^5 



2.5x10 



36 



Deharveng et al " ff§88 l 



Blair LonI| (|l997|>, opt, identified 



Deharveng et al. ( 1988 1 



Pannuti et al. 



(2000) 



Pannuti et al. 



(2000) 



air &; 



Deharveng et al, 



(1988) 



Blair & Long ( 1997|), opt, identified 



"[beharveng et al. ( 1988 ) 



optical sources are labeled. 

Two of the X-ray sources (4 and 8) were coincident 
with a saturated foreground star and the disk of a back- 
ground galaxy (with one visible point source nearly coin- 
cident with the X-ray position), respectively. Since both 
of these optical sources failed to meet the HST quality 
criteria of a stellar object in NGC 300, they were not in- 
cluded in the HST star catalogs. We used the IRAF task 
imexamine to estimate the magnitudes of the foreground 
star, the whole background galaxy, and the point source 
within the X-ray error circle centered on the disk of the 
background galaxy. The foreground star and background 
galaxy magnitudes were found to be consistent with those 
available in SIMBAD. Multi-band photometric measure- 
ments for each of the optical counterpart candidates to 
the X-ray sources are gi ven in Table [16] Source 4 is dis- 
cussed further in ^4.3.11 



4.3. X-ray /Optical Source Classifications 

The detection of optical counterpart candidates for 
each X-ray source with overlapping HST fields allows 
some sources to be assigned a source classification based 
on their multiwavelength properties. However, many 
X-ray sources are faint or do not have optical cover- 
age, making their physical nature ambiguous. We there- 
fore attempt to combine statistical quantities (i.e., using 
the radial source distribution, location in spiral arms, 
inter- arm regions, or background regions) with individ- 
ual source properties (variability. X-ray spectral shape, 
and optical colors) to determine the most likely source 
classification. 

We first consider the statistical properties of each X- 
ray source. For example, an X-ray source with a galac- 
tocentric radius less than 4 kpc that is coincident with a 
spiral arm is likely to be associated with NGC 300, while 
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Sources 32, 34 




Source 77 



Source 81 



Source 85 




Fig. 6. — Ho; emission from each of our SNR candidates (except source 46, for which Ho; imaging is not available). Each postage stamp 
is 0(5x0(5, centered on the apparent Ha emission. Crosses indicate the Chandra X-ray source position (determined by the 0.35-8 keV 
emission). The green contours show the 0.35-1 keV X-ray emission. 



an X-ray source at a large galactocentric radius that is 
not obviously associated with NGC 300 is likel y to be 
a background object. As described in Section |3.4[ we 
expect the majority of X-ray sources at radii beyond 4 
kpc are not associated with the stellar disk of NGC 300. 
We then examined ground-based, optical i?-band and Ha 
images to identify spiral arms and inter-arm regions in 
NGC 300, and define background regions not obviously 
associated with NGC 300. All sources with r > 4.5 kpc 
{0 > I'.l) were located in 'background' regions accord- 
ing to the ground-based imaging, making these sources 
probable background AGN. 

We next consider the individual X-ray properties of 
each source. Long-term X-ray variability is observed for 
both AGN and XRBs. However, SNRs and some HMXB 
systems do not sh ow significant X-ray variability (i.e., 
NGC 300 X-1; see [Binder erar||2011al and references 



therein). We therefore do not discount the possibility 
of an XRB nature for no n- variable sources. The data 
on galactocentric radius, location relative to NGC 300, 
and X-ray variability is then combined with the prelimi- 
nary source classificatio ns de fined using our hardness ra- 
tio analysis (see Section 3.2). For sources without over- 



lapping HST exposures, we then assign a likely source 
classification based on our X-ray observations. For ex- 
ample: a source with r = 1 kpc, located within a spiral 
arm, no significant long-term X-ray variability, and a soft 
X-ray spectral shape is classified as a likely SNR. 

To determine a likely source classification for X- 
ray sources with overlapping HST fields, we generated 
color-magnitude diagrams (CMDs) of each HST field 
and candidate optical counterparts were classified as 
falling either on the main sequence (MS) or red giant 
branch (RGB). Additionally, we constructed X-ray color- 




Fig. 7. — Optical RGB-rendered finding charts at the location of each X-ray source, indicated by a cross. White circles show the radius 
of the 90% Chandra error circle for each source. Candidate optical counterparts are indicated with red lines. X-ray sources with multiple 
candidate optical counterparts have their optical counterparts labeled 'a,' 'b,' 'c,' etc. 
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TABLE 16 

Optical Counterparts to Chandra X-ray Sources 



Source 
Name 
(1) 


(arcsec) 
(2) 


F435VF 
(3) 


F475Ty 

(4) 


F555W 
(5) 


F606W 
(6) 


F8UW 
(7) 


fb 

(10-1^ erg cm-2) 
(8) 


log(fx/fvY 
(9) 


4 (galaxy) 
4 (point source) 


1(^46 
0^18 


21.9 
24.7 




20.3 
23.8 




18.5 
23.7 


12.0 


-0.24 
1.17 


5a 
5b 
5c 
5d 


0^22 
0^/27 
0^/62 
0^70 


25.45d=0.05 


26.99±0.10 
25.51±0.03 
26.86±0.09 
26.24di0.05 


25.47±0.05 


26.39±0.06 
25.60±0.03 
26.17±0.05 
26.14±0.05 


25.47±0.08 
25.53±0.07, 25.86±0.08 
25.42±0.05 
26.52di0.18 


2.1 


1.44 
1.07 (1.12) 
1.35 
1.34 


7a 
7b 
7c 


0^/39 
0^42 
0('55 


26.91ib0.14 
27.13±0.18 
24.04ib0.02 




26.65±0.12 
26.93dz0.17 
24.13±0.06 




26.55di0.14 
26.13ib0.10 
24.25ib0.03 


1.6 


1.42 
1.53 
0.41 


8 


0'/04 


19.7 




18.5 




16.6 


3.9 


-1.44 


9a 
9b 


0'/34 
0^/39 


27.02±0.15 
25.98±0.07 


26.27±0.05 
25.43zt0.03 


25.48±0.06 
24.75±0.03 


25.09±0.02 
24.40±0.01 


23.98±0.04, 24.04±0.02 
23.51ib0.04, 23.53±0.02 


2.7 


1.19 (1.04) 
0.90 (0.76) 



Note. — "The positional offset from the Chandra X-ray source location. 
^X-ray flux is unabsorbed in the 0.35-8 keV energy band. 

^fy is calculated from either the FbbW or F6061^ magnitude, depending on which is available for a given observation. 
Only the flrst flve X-ray sources with optical counterpart candidates are listed. 
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Fig. 7. — (continued) Optical RGB-rendered finding charts at the location of each X-ray source, indicated by a cross. White circles show 
the radius of the 90% Chandra error circle for each source. Candidate optical counterparts are indicated with red lines. X-ray sources with 
multiple candidate optical counterparts have their optical counterparts labeled 'a,' 'b,' 'c,' etc. 



magnitude diagrams (XCMDs), where the optical color 
was plotted as a function of log{fx/fv)' Many authors 
have utilized optical and X-ray colors in discri minat- 



Horn- 



in^ between XRBs a nd back fflound AGN (e. 
schemeier et al.|2QQll [S hty: 

tably, 'iVLcGo wan et al.| ( |2U0 8) found that known pulsars 
and HMXBs located in the Small Magellanic Cloud oc- 
cupied a region of the XCMD bounded hy B — V < 
and fx/fv < 1- 

Optical counterpart candidates were then compared 
to the XCMD generated for each HST field. Sources 
with blue optical colors {F435W - F555W or F47bW - 
F606W < 0.2) and log{fx/fv) < 1 are likely to match 
a young massive main sequence star and are therefore 
designated as HMXB candidates. Sources falling in the 
other region of the XCMD are associated with either 
red giant branch (RGB) stars or background galaxies, 
and are therefore classified as either AGNs or LMXBs. 
For illustration, Figure [8] presents a sample CMD and 
XCMD for the field 'NGC300-5' (top panel) as weh as a 
CMD and XCMD of our X-ray/optical source candidates. 
The sample CMD was generated using the HST FASbW 
and FS14W catalog magnitudes for all optical sources in 
the 'NGC300-5' field. For illustrative purposes, we have 



created a sample XCMD using the mean X-ray fiux of 
sources in the NGC 300 catalog (~ 5 x 10"^^ erg s"^ 
cm~^, corresponding to 2 x 10^^ erg s~^ at the dis- 
tance of NGC 300). Using this X-ray flux, we calcu- 
late a log{fx/fv) for every optical source in the HST 
'NGC300-5' field. The bottom panels show the magni- 
tudes, colors, and X-ray/optical fiux ratios for the can- 
didate optical counterparts only. Regions corresponding 
to likely main sequence (MS) and RGB stars are labeled 
on the CMD, and the regions defining likely HMXBs and 
AGNs are labeled on the XCMD. 

In order to combine all the relevant X-ray and opti- 
cal information available from our Chandra and HST 
observations, we constructed an IDL routine called 
xclass.pro to automatically process each X-ray source. 
The routine reads in an input source file containing all 
X-ray and optical information (if available), and returns 
a preliminary source classification. The routine addition- 
ally tracks which classification channels were utilized by 
the input sources. The xclass.pro routine is publicly 
availabl^D 

^ Se e |http : //www . astro . Washington . edu/users/bbinder/l 
|xclass"7] 
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Fig. 8. — Top row: a sample color-magnitude diagram (CMD; left) and X-ray color-magnitude diagram (XCMD; right). The MS and RGB 
regions are marked on the CMD, and the HMXB and AGN regions are marked on the XCMD. Bottom row: the locations of our candidate 
optical counterparts as either MS or RGB stars on a CMD (left), and an XCMD (right) showing the actual values of the X-ray/optical flux 
ratio for our candidate optical counterparts to each X-ray source. 

in our X-ray source catalog to be ~50. The similari- 
ties in both the X-ray and optical properties of AGN 
and LMXBs (especially at the distance of NGC 300) is 
likely to cause confusion in separating these two popula- 
tions. With 47 AGN and 11 LMXB candidates obtained 
using our source classification scheme, we estimate our 
X-ray /optical classified \ogN-\ogS distribution to differ 
from the statistical by less than 10%. Our source clas- 
sification scheme correctly classifies both NGC 300 X-1 
(source 80) and SN 2010da (source 73) as HMXBs, which 
were previously identifi ed as HMXBs in earlier works (see 



The source classification reported by the xclass rou- 
tine was then compared to the classification determined 
by our hardness ratio analysis. In some cases, the multi- 
wavelength and hardness ratio classifications disagreed. 
This was most commonly seen with sources classified as 
XRBs by our hardness ratio analysis but with multiwave- 
length properties consistent with a background AGN. In 
these cases, we favor the AGN classification. Sources 
that are initially classified as AGN but which exhibit 
a soft X-ray spectrum (i.e., classified as 'SOFT' in our 
HR analysis) may be potential LMXBs. Only one of 
our HR-classified SNRs shows X-ray and optical char- 
acteristics inconsistent with a SNR origin, and may be 
an unusually soft background AGN or a foreground star. 
Our final source classification yields 46 background AGN 
candidates, 26 HMXB candidates, 11 LMXB candidates, 
10 SNRs, and one foreground star. A summary of all 
X-ray and optical data and our final source classification 
is present in Table pTl 

We use the [CappeTluti et al.| ( |2009| ) AGN logA/'-logS' 
distribution to estimate the expected number of AGN 



Binder et al.||2011b|a[ ). We are therefore confident that 
our source classitication scheme provides the most reli- 
able indicator of the physical nature of the NGC 300 
X-ray sources possible with the available data. 

4.3.1. Source 4' Foreground or Background? 

As can be seen in Figure [7| X-ray source 4 is coincident 
with the disk of a background galaxy. This source was 
not detected by XMM-Newton, suggesting a large bright- 
ening by a factor of > 22. The source has hardness ratios 
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TABLE 17 
X-RAY Point Source Classifications 



Source 


J.CL 


Location^ 


Variability^ 


In HST 


# Optical 


Location 


Location 


Prelim. 


ID. from 


Final 


No. 


(kpc) 


in NGC 300 




footprint? 


Counterparts 


on CMD 


on XCMD 


ID 


HRs 


ID 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1 


5.2 


bkg? 


none 


no 








AGN? 


ABS 


AGN 


2 


6.7 


bkg? 


long 


no 








AGN 


XRB 


AGN 


3 


5.5 


bkg? 


none 


yes 


2 


RGB 


AGN 


AGN 


SOFT 


LMXB? 


4 


6.7 


bkg? 


long 


yes 


1 






Galaxy 


ABS 


Galaxy/FS? 


5 


4.6 


bkg 


none 


yes 


10 


MS, RGB 


HMXB, AGN 


HMXB? 


XRB 


HMXB 



Note. — "Inclination-corrected galactocentric radius; the X-ray point source population at r > 4 kpc is expected to be dominated by 
background AGN. 

^Location in NGC 300 was determined by looking at ground-based i?-band and Ho; imaging. Locations are described as spiral arm, inter- 
arm regions, or background. X-ray sources that fall beyond the extent of the optical images are labeled as 'background?' 
^Variability is classifi ed as 'lonp: '-terra , 'rapid,' or 'none.' 

B. 

^NGC 300 X-1, see Bmder et al. {2GTldi,\. 



^SN 2010da, see Binder 



as lon g 
et al.ljs 



201IW 



Only the first five entries are shown.' 

consistent with those of other XRBs or background AGN. 
We consider three possible scenarios for the physical na- 
ture of source 4: (1) a background, high-redshift quasar, 
(2) an ultra-luminous X-ray (ULX) source located within 
the disk of the background galaxy, and (3) a flaring fore- 
ground M star. 

We flrst consider the high-redshift quasar interpreta- 
tion. Source 4 would have a typical quasar X-ray lumi- 
nosity of ~ 10^^ erg s~^ at a luminosity distance of ~ 26 
Gpc, roughly 2: ~ 3. However, high-redshift quasars typ- 
ically exhibit variations in X-ray luminosity on the order 
of ~ 15%, and extreme changes (such as seen in source 
4) are rare (Gibson & Brandt 2012). 

To evaluate the likelihood of this source being a ULX, 
we assume the background spiral galaxy has a physical 
diameter of ~22 kpc. We measure the angular diameter 
of the galaxy to be ^6^5 across, implying the distance 
to the host galaxy is ~700 Mpc. At this distance, source 
4 would have a 0.35-8 keV luminosity of Lx ^ 7 x 10^^ 
erg s~^. For a source radiating at the Eddington lumi- 
nosity, this implies a mass of ~5000 Mq. The factor of 
~ 22 brightening observed in this source is not typical of 
ULXs, with the notable exception of HLX-1 in the galaxy 
ESO 243-49 which has exhibited variations in X-r ay lu- 
minosity of a factor of ~20-50 (Lasota et aL]|20fT ). The 
high luminosity of HLX-1 is interpreted as mass accre- 
tion via Roche-lobe overflow onto a ~ 10^ Mq BH. We 
therefore flnd the ULX interpretation plausible. 

Young, active M dwarfs in the Milky Way undergo X- 
ray flares, and may provide a natural explanation for 
the factor of ^ 22 increase in X-ray luminosity exhib- 
ited by source 4. Assuming a typical M-dwarf My ^ 10 
and our estimated my ~ 23.8 for source 4 implies a 
distance of ~ 60 pc, yielding an X-ray luminosity of 
~ 5 X 10^^ erg s~^. The implied X-ray luminosity and 
color (F435W-F555W = 0.9) is typical of a flaring late- 
type star (ty pically 10^'^"^^ erg s~^ with B — — 
Singh et al. 1999). Using XSPEC, we simulate the spec- 



trum of a stellar corona using the MEKAL model with sim- 
ilar signal-to-noise as source 4 and estimate the hardness 
ratios produced. We flnd the hardness ratios of source 
4 to be consistent with the predicted hardness ratios, 
although at low signal-to-noise the spectrum cannot be 
distinguished from a power law. While we cannot rule 
out the possibility of an AGN or ULX origin, our obser- 



vations are most consistent with being a flaring, late-type 
star within the Milky Way. 

4.4. Comparison to Other X-ray Point Source 
Populations 

It is interesting to compare the properties of the 
NGC 300 point source population to other galaxies with 
a range of SFHs, morphologies, and stellar masses (M^). 
Since HMXBs serve as a tracer of young stellar popu- 
lations, their numbers sh ould correlate with the recent 
SFH of their host galaxy ( Grimm et al. [2003 ) , while the 
number of longer-lived LMXBs should correspondingly 
correlate with t he total stellar mass of the host galaxy 
( |Gilfanov||2004| ). Typically, the ratio L2-iokev/SFR is 
taken as a measure of the HMXBs in a galaxy, while 
LMXBs are measured through L2-iokev/^*- Both these 
ratios are found to correlate with th e morphology of the 
host galaxy (Grimm et al. 2002} [Gilfanov 2004). 

We estimate the 2-10 keV luminosity for our X- 
r ay /optically classifled HMXB candidates to be ^ 7 x 
]^q37 g-i g^j^^^ using equation (22) of Grimm et al. 



(|2003j), we derive a predicted SFR for NGC 300 of 0.12 
Mq yr-^ This predicted SFR is in excellent agreement 
with the average re cent SFR of 0.12±0.05 Mq yr"-*^ mea- 
sured by Gogarten et al. (2010). Normalized to the SFR, 
the HMXBs in NGC 300 emit - 6 x 10^^ erg ^'^/{Mq 
yr~^) in the 2-10 keV band, c onsis tent with the v alue 
measured for the Milky Way by 'Grimm et aL| (]2002|) . 

We flnd luminosity of our LMXB candidates normal- 
ized to the stellar mass of NGC 300 to be ~ 9 x 10^'^ erg 
s~^ Mq a factor of ^5 lower than is observed in the 
Milky Way (IGrimm et al.""2002). This value of LMXB 
luminosity normalized by stellar mass is a lower limit, 
as our source classiflcation scheme is unable to discrim- 
inate very well between LMXBs and background AGN. 
We can provide an upper limit to LMXB luminosity per 
unit stellar mass by assuming all AGN candidates are 
actually LMXBs - in this case, the luminosity normal- 
ized by the stehar mass of NGC 300 is ^ 6.5 x 10^^ erg 
s~^ M0 NGC 300 lacks the prominent nuclear bulge 
present in the Milky Way (where older LMXBs dominate 
the X-ray source population), and the gas in NGC 300 
(from which new HMXBs could be for med) accounts for 
^26% of the baryonic mass of NGC 300 (Westmeier e t al.| 
2011). Only ^ 15% of the baryonic mass of the Milky 
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Way ( Flynn et al.|2006 ) consists of gas that could poten- 
tially lorar^tarsTv^Tven the differences in stellar mass, 
gas content, morphology, and star formation history of 
NGC 300 and the Milky Way, this observed differences 
in LMXB populations is not surprising. 

When we consider the total XRB population, we find 
the number of X-ray sources above 10^^ erg s~^ i^x) 
and the total X-ray luminosity {Lx^ in units of 10^^ 
erg s~^) normalized by the stellar mass (in units of 
10^1 Mq) of NGC 300 to be Nx/ M^, = 93 dz 10 and 
Lx/M^ = OAT ± 0.05, respectively. |Gilfanov| ( |2004[ ) re- 
ported similar ratios of Lx/M^ for other late-type galax- 
ies, although the average value of Nx/M^ for the same 
galaxies was found to be ~ 160 ±50. It should b e noted 
that the comparison between the Gilfanov'(2004) sample 
and NGC 300 comes with two significant caveats: first, 
th e limiting X-ray luminosities of the galaxies examined 
by |Gilfanov| ( |2004| ) were > 10^'^ erg s"\ and the total 
number of X-ray sources above 10^^ erg s~^ often had to 
be estimated; second, the stellar masses of the galaxies 
sampled were significantly larger than NGC 300 (often 
by more than a n order of magnitude). Therefore, the 
ratios derived by Gilfanov (2004") may require modifica- 



tions for deeper X-ray data (below 10^^ erg s ^), and 
for galaxies with stellar masses comparable to NGC 300. 

5. THE logA^-logS' RELATION AND THE X-RAY 
LUMINOSITY FUNCTION 

Numerous studies have shown a correlation between 
the shape of the XLF of a po int source population and 
the SFH of the host galaxy (|Grimm et aL]|2003| lEra 



cleous et al.|[2Q06p . Those galaxies with a high degree ot 



recent ISl^' have their X-ray emission dominated by young, 
luminous HMXBs, while older stellar populations show 
LMXBs as the dominant source of X-ray emission, al- 
though this relationsh ip may not be mon otonic over a 
few hundred Myr (Belczynski et al. 2004). Using the 
Chandra X-ray point source catalog and overlapping HST 
exposures, we aim to directly test the XLF slope-age cor- 
relation in NGC 300. 

The cumulative number of sources above a given flux 
(A^(> 6'), where S is in units of erg s~^ cm~^) found in a 
survey covering a total geometric area A can be evaluated 
as: 



N{>S)= 



1 



i,Si>S 



A{Si) 



(4) 



A^(> S) has units of sources deg~^ and is weighted by 
the survey area over which each source could have been 
detected. The A{S) factor would account for survey 
completeness if the sources were distributed uniformly, 
although we do not expect this to be true for sources as- 
sociated with NGC 300. We terminate the low flux end of 
the distribution at the 90% limiting luminosity for each 
energy ba nd c onsidered. The sensitivity maps calculated 
in Section 2.4 allow the area function A(yS) (the area over 
which a source with flux S would be detectable) to be 
directly evaluated for each source. 

We use two approaches to construct the XLF for the 
NGC 300 point source population. First, we calculate 
the cumulative logA^-logS* distributions for the soft (0.5-2 
keV) and hard (2-8 keV) bands using all 95 X-ray sources 



detected in our Chandra o bservation. Using the b ack- 
ground AGN distribution of Cappelluti et al. (2009), we 
then statistically correct for the anticipated background 
contamination - we hereafter refer to these XLFs as our 
"statistical XLFs." 

The statistical XLFs are shown in Figure [9] For each 
energy band, we model the \ogN-\ogS distribution as a 
power law. For the 0.5-2 keV distribution, we find a best- 
fit power law index of 7 = 2.03 ± 0.10; there is a small 
bend in the XLF at ^ 1.6 x 10^^ erg s"\ although it 
is not statistically significant. Our 2-8 keV \ogN-\ogS 
distribution power law fit yields an index of 7 = 1.91 ± 
0.05. 

We next compute the cumulative \ogN-\ogS distribu- 
tion for objects classified as XRBs from our X-ray and 
optical classification scheme, hereafter referred to as our 
"X-ray /optically classified XLFs." We find the \ogN- 
\ogS distributions of our X-ray /optically classified XLFs 
agree with the statistical XLFs, especially in the 2-8 
keV energy range (due to less contamination by soft X- 
ray emitters such as SNRs). Figure 10 shows the X- 
ray /optically classified XLFs (red) ana our statistical 
XLFs (black) in both soft and hard bands. Due to their 
similarity in X-ray spectral shape and optical colors, the 
mis-classifications are likely to be LMXBs that are in- 
correctly identified as background AGN (or vice versa), 
while our HMXB sample is likely to have significantly 
fewer mis-classifications. The XLFs of X-ray /optically 
classified AGN are also consistent with the |Cappelluti| 
et al. (2009) AGN distribution. 

We tit the X-ray /optically classified XLF in both the 
soft and hard bands with power laws. The 0.5-2 keV 
XLF has a slope 7 = 2.42 ± 0.12, while the 2-8 keV XLF 
has 7 = 1.80 ± 0.07. A low-significance kink is observed 
in the soft XLF, similar to the feature observed in the 
soft statistical XLF. Our fits are consistent with those 
derived from the XMM-Newton point so urce catalog of 
NGC 300 above 1.5 x 10^^ erg s"^ (see|Carpano et al.' 



[2005^, their Fi gure 5). All the XLF slopes derived here 
are typica l of spiral galaxies with recent star formation 



activity ( [Fabbiano 2006 Colbert et al. 2004 Grimm et al. 
2003). 

To examine whether the bend in the soft XLF near 
3xl0~^^ erg s~^ cm~^ (equivalent to ~ 1.6 x 10^^ erg s~^ 
at the distance of NGC 300) may be due to a specific X- 
ray population we constructed logA^-logS* distributions 
for HMXBs, LMXBs, and SNRs separately. Figure [ll] 
shows the soft logTV-logS' distributions for each of these 
three groups, as well as the statistical XLF. We find that 
the HMXBs (black line) are the only population without 
a significant number of sources at fainter luminosities. It 
is possible that there exists a substantial population of 
lower luminosity, wind-fed HMXBs that would require 
deeper observations to detect at high significance. Al- 
ternatively, the change in the XLF slope may be due 
to XRTs that are observable only during their more lu- 
minous outbursting phase (and fall below our detection 



threshold during quiescence). Table 18 



summarizes our 



best-fit power law slopes for each of our XLFs. 

5.1. Radial logN-logS Distributions 

We next searched for radial trends in our NGC 300 X- 
ray source population. We divided our X-ray /optically 
classified XRB candidates into two radial bins: an "in- 
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Fig. 9. The cumulative 0.5-2 keV (left) and 2-8 keV (right) X-ray logA^-logS' distributions for the N GC 300 X-ray sources. The 
dot-dashed lines indicate the 90% limiting flux for our survey in each band. The AGN logA^-logS* model from |Cappelluti et al.| ( [2b09| > was 
assumed (red). The AGN-subtracted logA^-logS* distribution is shown in blue. 
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Fig. 10. — The 0.5-2 keV (left) and 2-8 keV (right) \ogN-\ogS distributions for our background-subtracted total \ogN-\ogS distributions 
(blue) and our X-ray/optically classifled XRB candidates (cyan). 



TABLE 18 

Power Law Fits to the NGC 300 logA^-logS' Distributions 



logA^-logS* Distribution 


0.5-2 keV 


2-8 keV 


Statistically XLF 


2.03±0.10 


1.91±0.05 


X-ray /optically Classifled XLF 


2.42±0.12 


1.80±0.07 


X-ray/optical XLF, r < 2.7 kpc 


1.68±0.31 


1.63±0.18 


X-ray /optical XLF, 2.7 < r < 5.4 kpc 


2.65±0.55 


1.23±0.07 


HMXBs 


1.86±0.19 


1.56±0.09 


LMXBs 


2.73±0.30 


1.69±0.29 


SNRs 


1.04±0.01 


N/A 



ner" region {d < 2.7 kpc) and "outer" region (2.7-5.4 
kpc) so that we could directly compare our logA^-logS* 
distributions to the radial SFHs measured by Gogarten 
et al.l (|2Q1Q[) Usi n^ the SFHs at different epochs (Gog- 
arten et al.||2UlU|), we assign an approximate, weighted 
average age to each radial bin. The radial \ogN-\ogS dis- 
tributions in both the soft and hard bands are shown in 
Figure 111 



We observe a preliminary correlation towards flatter 
2-8 keV XLF slope with increasing galactocentric radius 
and younger average age of the underlying stellar popu- 
lation. This trend is not observed in the soft band, which 
may suffer more severely from contamination from SNRs 
or foreground stars or from obscuration. We therefore 
consider the 2-8 keV XLF a more reliable diagnostic of 
the XRB population. In Figure [13) we compare our ra- 
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log Ln c 2 kev (el's s'^) at 2.0 Mpc 
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Fig. 11. The 0.5-2 keV logA^-logS' distribution for HMXBs (solid black), LMXBs (red), and SNRs (blue) as classified based on our 
X-ray and optical observations. The dashed line shows the statistical XRB logA^-logS* using the Cappelluti et al.^(2009 ) background AGN 
distribution. We find HMXBs are the only sub-group without a significant number of sources detected at tamter luminosities, potentially 
producing the break near 3 x 10"-*^^ erg s"-*^ cm"-^. 
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Fig. 12. — The 0.5-2 keV (left) and 2-8 keV (right) logA^-logS* distributions for X-ray/optically identified XRBs detected within the inner 
2.7 kpc (black) and outer 2.7-5.4 kpc (red) regions. 



dialXLF/age data to model predictions (Eracleous et al. 
2006 ). X-ray sources detected at smaller galactocentric 
radii (<2.7 kpc) show an average older age, ~83 Myr, 
and a steeper XLF slope (7 ~ 1.6) than those sources 
found in the outer regions of the NGC 300 disk (59 Myr 
in age with 7 ~ 1.2). 

Our X-ray /optical data of NGC 300 thus constitutes 
one of the first tests of such model predictions relating 
the nature of the X-ray sources to the underlying stellar 
population of a galaxy; additional XLF-SFH measure- 
ments will become available when the remaining four 
CLVS galaxies (NGC 404, NGC 55, NGC 2403, and 
NGC 4214) are analyzed in upcoming works. 



6. SUMMARY AND CONCLUSIONS 



We have presented the X-ray point source catalog for 
a new, deep Chandra observation of NGC 300 as part of 
the Chandra Local Volume Survey. A total of 95 sources 
were detected in NGC 300 down to a 90% limiting 0.35-8 
keV luminosity of 1.5 x 10^^ erg s~^. By comparing our 
X-ray source catalog to earlier observations of NGC 300 
performed by XMM-Newton^ we were able to place long- 
term variability constraints on each source and identify 
25 XRT candidates. We fit the radial X-ray source distri- 
bution to an exponential profile and were able to derive 
a scale length of 1.7±0.2 kpc, in good agreement with 
measurements taken at other wavelengths. We evaluated 
hardness ratios for each of our X-ray sources, and per- 
formed spectral fitting for those sources with more than 
50 counts. 

We utilized nine deep Hubble Space Telescope fields 



22 



Binder et al. 



Oh 

on 

X 
> 

M 

00 
I 




60 

Age (Myr) 



Fig. 13. — The slope of the cumulative 2-8 keV XLF as a function of stellar age. Black points indicate model predictions. The red cross 
indicates our data for radial bins less than 2.7 kpc, while the blue cross shows our data for X-ray sources from 2.7 - 5.4 kpc. Our combined 
Chandra/ H ST data of NGC 300 agree very well with predictions, and further analysis of the remaining CLV galaxies will provide further 
tests of these model predictions. 

to search for optical counterparts to each of our X-ray 
sources; 32 of our X-ray sources were observed in at least 
one of the overlapping HST fields. We developed a new 
source classification method which combined both X-ray 
and optical data to classify each source as a HMXB, 
LMXB, or background AGN, and have made our IDL 
routine which automates our source classification tech- 
nique publicly available. We classified 47 sources as back- 
ground AGN, 26 as HMXBs, 11 as LMXBs, 10 as SNRs, 
and one as a foreground star. 

The cumulative logA^-logS* distributions were pre- 
sented in both 0.5-2 and 2-8 keV bands and found to be 
consistent with those of other spiral galaxies with a domi- 
nant HMXB population. Additionally, we divided our X- 
ray sources into two radial bins and modeled the radially- 
resolved logT V-logS' distributions wi th power laws. Using 
the SFHs of Gogarten et al. (2010), we were able to as- 
sign an average to the underlying stellar populations, and 
we find the measured XLF slopes and stellar population 



ages to be in good agreement with model predictions. 
The remaining four CLVS galaxies (NGC 55, NGC 404, 
NGC 2403, and NGC 4214) wih provide additional con- 
straints on XLF-age models with upcoming works. 
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